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ABSTRACT 

We present abundance measurements from 26 emission-line selected galaxies at z ~ 0.6 — 0.7. By 
reaching stellar masses as low as 10* M©, these observations provide the first measurement of the 
intermediate redshift mass-metallicity (MZ) relation below 10^ M©. For the portion of our sample 
above M > l(f M© (8/26 galaxies), we find good agreement with previous measurements of the 
intermediate redshift MZ relation. Compared to the local relation, we measure an evolution that 
corresponds to a 0.12 dex decrease in oxygen abundances at intermediate redshifts. This result 
confirms the trend that metallicity evolution becomes more significant towards lower stellar masses, 
in keeping with a downsizing scenario where low mass galaxies evolve onto the local MZ relation at later 
cosmic times. We show that these galaxies follow the local fundamental metallicity relation, where 
objects with higher specific (mass-normalized) star formation rates (SFRs) have lower metallicities. 
Furthermore, we show that the galaxies in our sample lie on an extrapolation of the SFR-M, relation 
(the star- forming main sequence). Leveraging the MZ relation and star- forming main sequence (and 
combining our data with higher mass measurements from the literature), we test models that assume 
an equilibrium between mass inflow, outflow and star formation. We find that outflows are required 
to describe the data. By comparing different outflow prescriptions, we show that momentum driven 
winds can describe the MZ relation; however, this model under-predicts the amount of star formation 
in low mass galaxies. This disagreement may indicate that preventive feedback from gas-heating has 
been overestimated, or it may signify a more fundamental deviation from the equilibrium assumption. 



1. INTRODUCTION 

The balance between gaseous inflows, outflows, and 
star formation is a critical frontier in our understanding 
of galaxy evolution. Feedback caused by stellar winds, 
supernovae, and supermassive black holes is often used 
to explain a variety of observations, from luminosity and 
stellar mass functions, to the enrichment and reionization 
of the intergalactic medium. However, a c omplete physi- 
cal p i cture of these feedbac k processes (i.e. iMurrav et al.l 
120051: IHopkins et all I2012D is still debated. Continued 
efforts to provide new observational tests are essential. 

The correlation between galaxy stellar masses and gas- 
phase metallicities (the MZ relation) is one important 
probe of star formation feedback. The galactic out- 
flows that slow star formation, and the inflows that pro- 
mote it can also alter metallicities. On one hand, metal 
poor material, when accreted onto a galaxy, can lower 
the metallicity. On the other hand, supernova driven 
winds may remove metal enriched material from galax- 
ies. It has been known for many years that models 
which ignore inflows and outflows (i.e. closed-boxes), fail 
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to reproduce observe d abundance patterns in galaxies 
(jvan den Berghl [19621) . Indeed, in recent years, both 
analytical and numerical models have shown that out- 
flows (and sometimes gas accre tion) a re needed to ex- 
plain the MZ relation (Tremonti et al. | [2004: Dalcanton 
2007; Brooks et al. 2007; Finlator fc Davl [2 008: Erb 
^008; Dav e et al...2011b. .2012; .Peeples fc Shankar.,201lT 
[Dayal ct al.l 120131: IlTHv et al.ll2013D . However, to date, 
studies have not converged on the properties (rates, kine- 
matics, metal enrichment, halo mass dependence, and 
redshift evolution) of these gaseous flows. 

One avenue to better understand the physics that 
governs the MZ relation is to study low mass galax- 
ies (e .g. 'Lee et al.' '2006|; IZhao et al.r(2010l: IWuvts et all 
120121 : iBerg ct al. 2012). In these systems, galactic 
winds are especially effective at escaping the gravita- 
tional potential of their hosts, regulating star forma- 
tion and enabling enrichment of the intergala ctic medium 
(jOppenheimer et al.|[2009l: IKirbv et al.ll201lD . Hence, ob- 
servational constraints on low mass galaxies offer some 
of the most stringent tests of galaxy formation models. 

Outside the local universe, the MZ relation is 
poorly constrained at stellar masses below 10^ Mq 
(and 10"^° Mq for z > 1). Wh i le large spectro- 
scopic survevs (iLe Fevre et al.l 1 20041: iLillv et all 120071: 
iNewman et al.l[2012l ) have enabled abundance measure- 
ments of statistical samples out to z ^ 1, these sur- 
veys are typically limited to i? < 24. Hence, these 
intermediate-red shift MZ relations have been derived for 
M > l(f M© .(I.LiUv et al.l [20031: [Sa vaglio ct al. 2005] 



Lamare ille et al.l l2009t Zahid et al.l 1201 1; iCresci et al 









2012t [Moustak as et al.l[2C)llt ). Nevertheless, in the case 



of the Cosmic Origins Survey (COSMOS), the extensive 



Table 1 

DEIMOS FoUowup Observation Summary 



Date 



Mask Name 



Mask RA 
(J2000) 



Mask Dec 
(J2000) 



Mask PA 
(degrees) 



Slit PA 
(degrees) 



Slit Widths 



Exposure Time 
(hours) 



Seeing 

(") 



27 January 2011 D 10:00:22.97 02:09:28.9 85 90 1.5 

28 January 2011 F 10:01:11.46 02:10:27.4 106 90 1.5 
22 January 2012 M 10:00:24.25 02:05:19.8 85 90 1.5 

22/24 January 2012 L 10:00:22.56 02:15:29.8 13 1.2 

23/24 January 2012 Q 10:00:28.96 02:20:05.1 95 90 1.5 



6.5 
6.3 
4.9 
5.3 
6.8 



0.6 
0.8 
1.0 
0.9 
1.0 



Note. — Coordinates, position angles (PAs), exposure times, and seeing are given for each observed mask. O n masks D, F, 
M, and Q, the slit PAs a nd slit-widths were matched to the IMACS Venetian blind spectroscopic data (described in lDressler et all 
I2011al and IHenrv et all 2012). For mask L, the slits were rotated 90 degrees relative to the search data. 



broad and intermediate-band photometry allows reliable 
mass constraints an order of magnitude lower at inter- 
mediate redshifts. Therefore, spectroscopic foUowup of 
fainter galaxies can significantly extend the intermediate 
redshift MZ relation. 

In this paper we use an emission line selected sample 
to place new constraints on the low mass end of the MZ 
relation at z ~ 0.6 — 0.7. By drawing our sample pri- 
marily from the ultra-faint emission line objects that we 
have previously ide ntified wi th blind spe ctroscopy in the 
COSMOS fi eld (M artin et al . 2008; Dress ier et al.ll2011al : 
IHenrv et~ai.l 120121) , we obtain oxygen abundances for 
galaxies with stellar masses of 10* M© < M < 10^° Mq. 
In this manner, we provide the first constraints on the 
low mass evolution of the MZ relation, reaching stellar 
masses that are comparable to the limiting mass of local 
SDSS samples. 

This paper is organized as follows: in f}2] we describe 
our spectroscopic observations and the COSMOS imag- 
ing data that we use. In ij3]we describe our emission line 
measurements and stellar mass derivations. Then, in 12] 
we calculate the oxygen abundances of galaxies in our 
sample, discussing the various diagnostics that have been 
proposed to break the de generacy of the do uble- valued 
R23 metaUicity indicator (jPagel et al.lll979t ). In §5 and 
§6 we compare our MZ relation to previous derivations, 
and investigate the presence of a mass-metallicity-SFR 
fundamental plane. Finally, we compare to theoretical 
predictions of the MZ relation in §7. In this paper we 
use AB magnitudes, a Chabrier (2003) initial mass func- 
tion, and a ACDM cosmology with J^m = 0.3, f^A = 0.7 
and 7?o = 70 km s~^ Mpc~^. Throughout the text we re- 
port measurements of doublet lines: [O II] AA3727, 3729 
and [O III] AA4959, 5007. For the sake of brevity, we use 
the notation "[O II]" and "[O III]" to refer to both lines 
in the doublet, or, when appropriate, the sum of their 
fluxes 



2. OBSERVATIONS 



2.1. 



Target Selection & FoUowup Spectroscopy 

The emission line galaxies in the present sample were 
initially identified as part of our multislit narrowband 
spectroscop ic survey. T he obs ervations are p resented 
in detafl in 'Martin et alj (|2008D . [Dressier et all (|2011aD . 
and Henrv ct al. (201^. In brief, our design uses 
the Inamori- Magellan Areal Ca mera and Spectrograph 
(IMACS; iDressler et all l2011b[ ) on the 6.5m MageUan 
Baade Telescope at Las Campans Observatory. We used 
a Venetian blind slit mask, and narrowband filter cen- 
tered in the 8200A OH airglow free window. This method 



allows for the efficient selection of emission line galaxies, 
and reaches fluxes as low as 2.5x10^^* erg s~^ cm~^- 
a fact or of flve fainter than n arrowband imaging surveys 
(e. g. lKashikawa et al.ll201l"l ). 

Because of the faintness of the emission lines, and 
the narrow bandpass in our search data, we use fol- 
lowup spectroscopy with DEIMOS (Fa ber et all I2003D 
to identify the redshifts of the IMACS-detected galax- 
ies. The primary goal of these observations was to 
confirm Lya emitting galaxies and measure the faint- 
end slope of the L yct luminosity function at z = 5.7 
(jHenrv et al.ll201^ . However, for the foreground galax- 
ies at z ~ 0.6 — 0.7, these followup observations con- 
tain the [O III]AA4959, 5007, [O II]AA3727, 3729, and H/3 
emis sion lines that com prise the R23 metallicity indica- 
tor (jPagel et al.l 119791 ) . Three classes of these objects 
were included on followup slit-masks: 

1. Emission line galaxies for which continuum in the 
search data ruled out the Lya identification were 
chosen as MZ targets. To obtain the best possible 
spectrum, we searc hed the COSMOS photomet- 
ric redshift catalog (jllbert et al.l 120091 ) for nearby 
galaxies (< 2") with z ~ 0.6 — 0.7, and shifted 
the slit-positions to coincide with these matchetQ- 
In total, the final sample of 26 star forming galax- 
ies (described in Q contains 15 objects that meet 
these criteria. 

2. Additionally, galaxies were drawn from the [O III] 
-I- up narrowband exces s catalog which we derived 
in lDressler et al.l ()2011aD . The present star-forming 
sample contains 11 objects that were selected in 
this way. 

3. Finally, we also detect R23 from ultra- faint emis- 
sion lines which we originally considered Lya can- 
didates, but, upon the followup observations de- 
scribed below, we determined that the discovery 
line was [O III] A4959, [O III] A5007, or H/3. Six 
galaxies with measurable R23 fall under this clas- 
sification. 

In practice, the objects described in #3 are difficult to 
evaluate; their imaged counterparts are sometimes am- 
biguous due to their faintness and uncertain position 
within the blind-search slit. Therefore, in this paper we 

^ Because the COMSOS photometric redshifts include narrow- 
band photometry and use templates that include emission lines, 
the photometric redshifts are often very close to the spectroscopic 
redshifts for the emission-line selected galaxies in our survey. 
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Figure 1. The spectra of 34 galaxies witii complete coverage of the R23 metallicity indicator ([O III]+ [O II]/ H/3: IPagel et al.lll979l ) are 
shown. The spectra of 26 objects that comprise our MZ sample are shown in black, and are sorted by stellar mass with higher masses at 
the top (the same as their appearance in Tables [2] and [Sj . The two AGN candidates are shown in grey at the top, and also continue the 
mass-ordered scheme. The Class 3 objects with ambiguous optical counterparts are shown at the bottom (in grey) in the same order as 
their appearance in the tables. All spectra are normalized by their flux in [OIII] A5007, so it is apparent that higher mass galaxies have 
relatively stronger H/3 and [O II] emission. The strong "absorption" trough near 4845 A in one object (980802) is a detector artifact. 



focus on the former two classes. For completeness, in 
Tables 2 and 3 we list the subset of these objects where 
R23 could be measured; however, we do not consider 
them further in this paper. 

The DEIMOS observations of these foUowup slit-masks 
were carried out in January 2011 and January 2012. A 
summary of observations for each followup mask is given 
in Table [TJ In four of the five masks, we used slit- widths 
and PAs that were matched to the search data (1.5" wide 
and 90° east of north). On mask L, in order to better 
locate objects detected through blind spectroscopy, we 
used a slit orientation which is orthogonal to the Venetian 
blind search slits. This method also allowed for narrower 
slits (1.2"). All observations were carried out using the 
the GG495 blocking fiher with 830G grating. Under this 
configuration, we achieved a spectral resolution of 3.7 A 
(2.9 A) for a source that fills the 1.5" (1.2") slit. The 
grating angle was chosen to give a central wavelength of 
7270 A, with blue coverage down to 5500A for most slits. 

The DEIMOS spectra were reduced using the DEEP2 
DEIMOS data reduction pipeline fCoopcr ct al. 2012), 
with an updated optical model for the 830G grating (P. 
Capak, private communication). The data were flux cal- 
ibrated using observations of several spectrophotometric 
standard stars, taken through a 1.5" slit at the parallac- 
tic angle. The stars use d were G191B2B, GPSO, Feige 
66, F eige 67, and Hz 44 (jMassev fc Gronwalll[l990l : \OM 
Il990t ). and the data for each was taken from the ESO 
spectrophotometric standard star database. Sensitivity 
functions derived from these stars differ systematically 
between observations, with offsets up to 30%. After shift- 
ing the lower sensitivity functions to match the highest 
one, the observations agree at the 2-3% level, indicating 
an excellent relative calibration. 

Finally, we have verified that the effects of differential 
atmospheric refraction have a negligible impact on our 
measured line flux ratios. In order to quantify possible 
differential slit losses, we calculate the component of the 
atmospheric refraction that falls perpendicular to the slit 
in each observation (frame) of each mask. The worst case 
occurs on mask L, where the slits are narrowest and the 
slit-PA was the furthest from parallactic. For these data, 
the refraction (perpendicular to the slit) between 6000 A 
(near [O II]) and 8000 A (near [0 III] and H/3) ranges 
from 0.04 to 0"33 (|Fillipenkolll"98a ). Therefore, since the 
guiding was done in the R-band (in between the observed 
wavelengths of the emission lines), we estimate the slit- 
losses on the red and blue parts of our spectra differ by 
no more than 4% (for 1" seeing). In most frames, the 
effect is even smaller, so we conclude that no systematic 
correction is needed to interpret our emission line ratios. 

2.2. COSMOS Imaging Data 

In order to derive stellar masses (SJ3]), we use the 
wide range of imaging data provided by the COSMOS 
team. A more detailed descr iptio n of the catalogs ca n be 
found in lCapaketHI (|2007| ) and lllbert eFaD (|2009| ). In 
summary, the photometric data include GALEX imag- 
ing, a wealth of ground-based, broad-band optical and 
near-infrared imaging, and the four Spitzer/IRAC bands. 
Furthermore, our SED fits are improved by the inclu- 
sion of intermediate-band imaging, and updated Sub- 
aru/SuprimeCam z'-hand data (from observations made 



with the new, fully depleted, red-eflicient CCDs). Since 
the optical and near-infrared photometry are performed 
in 3" apertures on data with homogenized PSFs (1.5" 
FWHM), we apply a point source aperture correction of 
-0.28 magnitudes to all optical and near-infrared bands. 
For the IRAC data we use 3.8" diameter apertures, and 
apply point source aperture corrections of -0.29, -0.33, 
-0.51, and -0.59 magnitudes to bands one through four. 
No aperture corrections are applied to the GALEX data, 
as these magnitudes were derived from fits to the PSF. 
Finally, we note that th e recommended zero point off- 
sets have been applied (jCapak et al.l l2007t lllbert et"all 
l2009f ). and a Gal actic foreground exti nction correction 
is made using the iCardelli et al.l (|1989[ ) extinction curve 
with E(B-V) = 0.019. 

3. ANALYSIS 
3.1. Emission Line Measurements 

In total, we identify 47 galaxies at z ~ 0.6 — 0.7; of 
these 34 h ave complete cov erage of the R23 ([O III]-|- 
[O 11]/ H^g jPagelet al.ll979D metalhcity indicator. Their 
spectra are shown in Figure [TJ The remaining 13 galaxies 
at these redshifts either have essential lines lost in the 
OH airglow spectrum, or are class #3 objects (described 
above), where the observed spectrum is especially faint 
because the galaxy may not be centered in the slit. Only 
six class #3 objects (without mass measurements) have 
R23 measurements; their spectra are shown in grey at 
the bottom of Figure[T] Excluding these objects leaves 28 
galaxies with both mass and R23 constraints. In §3.3 we 
show that two of these galaxies may contain active nuclei, 
so we adopt a final sample that includes the remaining 
26 star-fornfing objects. 

Emission line fiuxes are measured for all of the z ^ 
0.6 — 0.7 galaxies that we identified. Gaussian profiles are 
fit to each line, and uncertainties are determined through 
a Monte Carlo simulation where noise is added to the 
spectrum and the fits are repeated. For [O II], the two 
doublet components are fit simultaneously. As an exam- 
ple. Figure [2] shows a typical spectrum and the Gaussian 
fits to the emission lines. 

To facilitate a correction for H/3 stellar absorption, we 
also measure emission line equivalent widths. Because it 
can be difficult to assess the uncertainty on the contin- 
uum flux density when it is measured spectroscopically, 
we compare two methods. On one hand, we measure the 
equivalent widths directly from the spectra. In 25/28 ob- 
jects in our sample, we detect continuum, although it is 
sometimes weak. As a second method, we compare the 
emission line fluxes from the DEIMOS data to the contin- 
uum under the emission line from the SED fit (described 
below). After accounting for a systematic offset (because 
the emission line fiuxes are subject to slit losses but the 
SED fits use total fluxes), the equivalent widths mea- 
sured by these different techniques agree to within 50%. 
We adopt this level of uncertainty on the H/3 equivalent 
widths. 

Finally, having measured the fluxes and equivalent 
widths of our emission lines, we apply a correction for 
the stellar absorption. Because the stellar absorption 
component is much broader than the emission line, this 
correction depends on the spectral resolution. At the res- 
olution of our data, the correction is approximately 1 A 
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in eq uivalent width (jCowie fc Barged 120081 : IZahid et alJ 
1201 11 ) . While the amount of absorption also depends 
slightly on the age of the stellar population, the main 
source of uncertainty in this correction is from the ob- 
served H/3 emission equivalent width, as discussed above. 
We propagate this error in our H/3 fluxes and the metal- 
licities that we infer in 341 



Figure 2. The DEIMOS spectrum of object 54.5+4-0.18 measures the R23 oxygen abundance indicator with high signal to noise. The 
Gaussian fits used to measure the emission line fluxes are shown in red. From left to right, the lines shown are [O II] AA3727, 3729, H/3, 
[O III] A4959, and [O III] A5007. This spectrum is representative of the fainter, lower mass galaxies in our sample (see Tables [2] and [5]| . 

subtracted from the r and i— band continuum flux den- 
sities. Emission from Ha, on the other hand, falls be- 
tween the z' and J bands, so it does not contaminate 
any of the photometry in the present analysis. Finally, 
we consider emission line contamination to the interme- 
diate band photometry. In these cases, the corrections 
will be larger and more uncertain. Therefore, we exclude 
the band at 624 nm that includes [O II] , and the band at 
827 nm that covers [O III] and H/3. The remaining ten 
intermediate bands should be relatively unaffected. 

The grid of stellar population parameters that we fit 
with FAST include: a set of exponentially declining 
star formation histories with e-folding times, r, ranging 
from 40 IVIyr to 10 Gyr; characteristic stellar population 
ages ranging from 50 JVIyr to the age of the u niverse at 
z ~ 0.6 - 0.7; and Ay = - 3 for a lCalzetti ef al. (200(1) 
extinction curve. Additionally, we use a lChabrierl (|200j) 
IJVIF with metaUicities of 0.004, 0.008, and 0.02 (solar). 
Supersolar metaUicities are excluded, because, as we will 
show in fj4l the galaxies in our sample mostly have sub- 
solar to solar gas-phase metaUicities. The derived stellar 
masses, SFRs, and visual extinction values are given in 
Table [3l Because of the intermediate band photome- 
try, the dust extinction constraints exclude much of the 
allowed parameter space. Therefore we use the SED- 
derived dust constraints to correct our emission line ra- 
tios, including these uncertainties in our error budget. In 
making this correction, we also account for the fact that 
nebular extincti on is 2.3 time s higher than stellar extinc- 
tion on average (jCalzetti et a l. 2000; Cresci et al. 2012). 
While the precise relation betwe en stellar and nebular 
extinction remains controversial (iCowie fc Barged [20081 : 
iCresci et all 120121: [W offord et a0l2013D, we adopt the 
"standard" iCalzetti et al.i (|2000il relation to allow com- 
parison with other studies. 

In ij6] and [7] we draw conclusions based on the SFRs 
that we have derived from these SED fits. In order to 
verify that the SFRs are reliable, we have carried out an 
assessment of the systematic uncertainties by comparing 
SFRs derived from [O II], H/3, and the SED fits. The re- 
sults, outlined in the Appendix, show that the agreement 
between the different diagnostics is good; systematic off- 
sets are smaller than 0.2 dex. This level of uncertainty 



3.2. Stellar Masses, SFRs, and dust constraints from 
SED fitting 

Stellar masses are determined by fitting template stel- 
lar populations to the broad- and intermediate-band 
photometry described in §2. For this task, we use 
FAST (Fitt i ng an d Assessment of Synthetic Templates; 
IKriek et al.l 120091 ). The choice of population synthesis 
templates can have an important impact on the proper- 
ties derived from SEDs, because the contribution from 
the (infrared-bright) thermally pulsing asymptotic giant 
branch (TP-AGB) stars differs from model to model. 
The most widely used templates are Bruzual & Chariot 
(2003; BC03). Maraston (2005) . and Chariot & Bruzual 
(2007; lBruzualll2007[ ): the latter two options provide a 
larger contribution from TP-AGB stars (and correspond- 
ingly smaller masses). Nevertheless, the proper contri- 
bution from TP-AGB stars remains a subject of debate 
(IKriek et all 120101 : iConrov fc Gunn 2010;' Zibctt i et al.l 
I2013f ). Therefore, in order to facilitate comparisons with 
the literature, we derive stellar masses using the BC03 
models. We have verified that the Chariot & Bruzual 
and JVIaraston et al. models produce stellar masses that 
are systematically smaller by approximately 0.1 dex. 

In order to assure accurate stellar population con- 
straints, it is also important to account for emission line 
contribution to our SEDs. (The stellar synthesis tem- 
plates that we fit to the SEDs do not contain nebular 
features. Failing to remove emission line contamination 
can result in in correctly inferred age s and stel lar masses; 
ISchaerer fc de"B arros 2009: Atc k etHI [20111 ) We use 
our line fluxes to calculate the small contribution from 
[O II] to the r-band, and [O III] and H/3 to the i-band. 
For the galaxies in our sample, we find that 1-5% of 
the r— band light and 2-15% of the i— band light can 
be attributed to emission lines. These contributions are 
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Figure 3. The [O III] / H/3 ratio as a function of stellar mass 
can be used to identify A GN, similar t o the fO III]/ H/3 vs . 
[N 11]/ Ha "BPT" diagram fHildwin, Phi llips, fc TerlevichI [TMil : 
IJun eau et al."201V). Red point show our emission line selected 
sample, while contours (with levels defined arbitrarily) show the 
SDSS data. The solid black lines show the demarcation between 
star- forming galaxies (below and to the left), objects with active 
nuclei (above and to the right), and composite objects (the region 
in between the black lines around 10^"- 10^^ I^©)- The two high- 
est mass objects in our sample, 55.5-3-0.60 and 73.5-1-5-0.54, fall 
outside of the z = star-forming locus, so we identify them as 
candidate AGN. 

does not affect our concfusions. 

3.3. Contamination from AGN 

fn order to measure tlie MZ relation of star-forming 
galaxies, it is important to ensure that the observed 
emission lines do not originate from an active galac- 
tic nucleus (AGN). The typical approach for low red- 
shift galaxies is to use the [O III]/ H/3 and [N 11]/ 
Ha emission l ine ratios (th e canonical BPT diagram; 
[Baldwin, Philli ps, fc Terlevi ch 1981). However, for the 
redshift of our sample, Ha and [N II] fall at infrared 
wavelengths. Therefore, in lieu of foUowup spectroscopy, 
we turn to an alt ernate diagnostic: th e Mass-Excitation 
(MEx) diagram ( Juneau et al.l [20Tll ). This approach, 
shown in Figure [31 uses stellar mass as a proxy for the 
Ha/[N II] A6583 ratio (relying on the IVIZ relation). The 
resulting diagram appears qualitatively similar to the 
traditional BPT diagram, with AGN falling towards the 
top and right of the plot. In this figure, the data for 
our galaxies (red points) are compared to that from the 
SDSS (contours), and the solid black line demarcates the 
difference between star-formi ng galaxies, AGN, a nd com- 
posite objects (as defined bv lJuneau et al.ll2011l ). While 
the evolution of IVIZ relation compromises somewhat the 
use of this method, we note that two of our galaxies are 
found in the area identified as possible AGN. Since metal- 
licity evolution works in the sense of shifting the proxy 
relation between Ha/[N II] A6583 and mass, it shifts the 
star-forming locus (and black threshold curves) to the 
right in Figure [3l Therefore, the two objects in the AGN 
or composite part of the sample may in fact be part of 
the non-AGN distribution; however, we take the conser- 
vative step of excluding them from further analysis. We 
list the measured and derived properties of the two can- 
didate AGN at the ends of Tables [2 and [3 

4. CALCULATING OXYGEN ABUNDANCES 

Interpreting metal abundance measurements requires 
that we account for systematic uncertainties. Dif- 



ferent methods for measuring strong oxygen abun- 
danc es yield results which are offset by up to 0.7 
dex CKewlev fc Elli sonll2008t iLopez-Sanchez etal.1 [20121: 
[Andrews & Martini 120131 ). On one hand, photoioniza- 
tion models that are used to derive theoretical cali- 
brations (i.e. Kewlev fc Dopital [2002) are often based 
on simplistic assumptions about H II region geome- 
tries and ionizing spectra (jvan Zee et al.l 120061 ). On 
the other hand, empirical methods, which correlate 
line ratios with H II region electron temperatures (i.e. 
iPettini fc P agcl 2004; Pilyugin fc Tliuan 2005), gener- 
ally give lower metal licities than, pho t qionization model- 
based calibrations (jMcGaughl 119911 : iKewlev fc Dopital 
I2002D . Because electron temperature measurements may 
be overestimated when temperature and density gradi- 
ents are present in H II regions, it is possible that the 
empi rical calibrations are biased towards low m etallici- 
ties (iPeimbert fc Coster"olll969: Stasinskal 120051 ). Alter- 
nativelv. iNicholls et al.l (|2012D have suggested that dis- 
crepancies between electron-temperature measurements 
and theoretical strong-line estimates can be explained if 
the electrons in H II regions deviate from an equilibrium 
Maxwcll-Boltzmann distribution. 

Ultimately, the differences between metallicity calibra- 
tions are still not fully understood, but the systematic 
offsets can be accounte d for using a set of tran sforma- 
tion equations given by IKewlev fc Ellison! (|2008D . In the 
sections that follow (except for §6), we use the calibra- 
tion from Kobulnicky fc Kewley (2004; hereafter KK04), 
which is the me a.n of two theore tical R23 calibrations 
in th e literature (|McGaughl 11991! and IKewlev fc Dopita! 
20021). In §6 we instead use the calibration from 



Maiolino et al.! (I2008D so that we may compare our data 



to results from iMannucci et al.! (|2011l ) 

4.1. High- Metallicity or Low? Determining the Branch 

of R23 

Metallicities derived from R23 can be degenerate, since 
this diagnostic is double valued. At high metallici- 
ties, the oxygen lines are weaker because cooling is ef- 
ficient and the H II regions have lower electron tem- 
peratures. On the other hand, at low metallicities, the 
overall decrease in oxygen relative to hydrogen also im- 
prints a decrease in R23 with decreasing metallicity. The 
"turnover" metallicity that demarcates the transition be- 
tween the upper and lower branches of R23 depends on 
the ionization parameter, and differs among the calibra- 
tions that are in the literature. For the typical ionization 
parameters of our galaxies, and the KK04 calibration, 
the turnover metallicity is around 12-|-log(0/H) ~ 8.4, 
which is reached when log(R23)^ 0.9 — 1.0. 

Many different methods have been suggested for break- 
ing the degeneracy between the upper and lower branch. 
The most rigorous of these is to compare to an alter- 
native metallicity diagnostic. IKewlev fc Ellison! (|2008f ) 
advocate for the use of the [N II] A6583/[0 II] or [N II] 
A6583 /Ha ratios to first provide a metallicity estimate 
that identifies the branch. However, this method is more 

^ The ionization parameter is defined as the ratio of the ionizing 
photon density to the hydrogen density. It can be written as U = 
Q/iirr^njjc, where Q is the ionizing photon rate, r is the radius 
of the H II region, n is the hydrogen density, and c is the speed of 
light. The parameterization g = [/ X c is also commonly found in 
the literature. 



challenging at z > 0.5 (requiring infrared spectroscopy 
to reach very faint [N II] AA6548, 6583 hnes), and impos- 
sible from the ground at z > 2. The R23 diagnostic, on 
the other hand, can be measu red from the ground out 
to z ^ 3 (jMaiolino et alj [20081 ). It is therefore impera- 
tive that we learn how to break the R23-metallicity de- 
generacy without observations of the [N II] AA6548, 6583 
lines. Here, we examine other methods for establishing 
the R23-branch that have been proposed in the litera- 
ture. 

In Figure m we investigate the branch of R23 by plot- 
ting four different diagnostics as a function of stellar 
mass. For comparison, we show the four H/3-detected 
galaxies (0.60 < z < 0.85) from the PEARS Survey 
(Probing Evolu tion and Reionization Spectroscopically; 
iXia et ani201^ . In the upper left panel, we show that 
R23 decreases with increasing stellar mass for most of 
our sample. Furthermore, three of the four galaxies from 
Xia et al. follow the trends from our emission line se- 
lected sample. This diagnostic paints a clear picture: if 
metallicity decreases as a function of decreasing mass, 
most of the galaxies in our sample must fall on the up- 
per branch of the R23 indicator. The two lowest mass 
galaxies shown in Figure |4] (1217842 from our sample, 
and 246 from Xia et al.) hint at a possible turnover at 
M r^ 8.0 — 8.5, indicating that these galaxies may fall on 
the lower branch of R23. 

Because it is not ideal to introduce a mass dependence 
on our metallicity derivations, we explore other emis- 
sion line diagnostics in Figure |4l First, we compare to 
the emission line ratios of [O III]/H/? and log([0 III]/ 
[O II]) = 032 as a function of stellar mass. These line 
rat ios have been propose d to break the R23 degeneracy 
bv iMaiolino et all (|2008t ). with low metallicity solutions 
being preferred when [O III]/H/3> 5 and 032 > 0.5. 
Indeed, Figure |4] shows that the galaxies above these 
thresholds also prefer to have higher values of R23 (and 
lower metallicities on the upper branch). Nevertheless, 
we still see very little evidence for a substantial popula- 
tion of galaxies extending to metallicities below the R23 
turnaround. The only object that passes bot h the 032 
and [O III]/H^ thresholds is object 246 from iXia et~all 
((20T2I) . 

The inference of uppe r branch metallic ities for 3/4 of 
the objects presented bv lXia et al.l ()2012|) contrasts with 
their lower branch assumption. These authors argued for 
lower metallicities because of the relatively high values of 
[O III]/H/3 and 032 for most of their sample. However, 
comparison to the R23 vs mass plot (upper left) shows 
that for at least three of these galaxies, higher metallic- 
ity solutions are more plausible. (Although the galaxies 
presented by Xia ct al. 2012 may prefer higher 032 and 
ionization parameters than the remainder of our sample, 
metallicity is not a strong function of ionization parame- 
ter on the upper branch, so their galaxies still follow our 
R23-stellar mass correlation.) 

Ultimately, caution is required when using [O III]/H/? 
and 032 to determine the branch of R23, because these 
quantities are dependent on the ionization parameter. 
Observations of high-redshift galaxies show that their 
line ratios are offs et to high values of [O III]/H/3 at 
a fixed [N II] /Ha (IShap lev et al.l[200a lErb et al.l [20061: 
iHainhne et al.|[2009 ). This offset is usually interpreted as 
evidence for a systematically high ionization parameter 



(jBrinchmann et all I2008D . suggesting that some galax- 
ies deviate from the local metallicity vs. ionization pa- 
rameter correlation. In summary, attempting to break 
the R23 degeneracy using quantities that depend on the 
ionization parameter could erroneously indicate lower 
branch solutions. 

In Figure[ll we also investigate the use of H/3 equivalent 
width to s elect low-metalli city galaxies. iKakazu et al.l 
(|2007D and iHu et al.l (|2009D have found that, for their 
emission line selected sample, galaxies are likely to have 
very low metallicities (indicated by detectable [O III] 
A4363 emission) when their H/3 equivalent widths are 
more than 30 A (rest). However, 8/26 galaxies in our 
sample meet this criterion, and we do not detect [O III] 
A4363 in any objects. What is more apparent in Figure 
[H is that a cut in H/3 equivalent width merely selects 
against galaxies with M > 10®° Mq. As we show in iJH 
the Hu et al. sample has, on average, higher H/3 equiv- 
alent widths and somewhat lower luminosities than the 
present sample of emission line objects. We infer that it 
is not straightforward to use H/3 equivalent width to dis- 
criminate between upper and lower branch metallicities. 

Guided by Figure [4| we conclude that, in the absence 
of other metallicity diagnostics, the correlation between 
stellar mass and R23 is the preferred method to break 
the R23-degeneracy. This exercise shows that for z ~ 
0.6 — 0.7 galaxies, the maximum value of R23 is reached 
between 10^° and 10*'^ M©. Therefore, we adopt upper 
branch solutions above M = 10®'^ Mq. The two lowest 
mass galaxies ( 1217842 from our sample and 246 from 
IXia et al.ll2012D are less certain. For 1217842 we adopt 
a metallicity in the turn around region (12 + log(0/H) 
= 8.4), with error bars denoting the upper bound of the 
high metallicity solution and the lower bound of the low 
metallicity solution. The lowest-mass object from Xia et 
al. (#246) is tentatively assigned to the lower branch 
because of its more extreme line ratios, although we cau- 
tion that the H/3 measurement could be in error because 
it is blended with [O III] in their low resolution data. 
FoUowup spectroscopy can better constrain the metallic- 
ities of galaxies near the R23 turnar ound by providing 
the [N II] /Ha and [N II] /[O II] ratios (jKewlev fc EUisonl 
l2008f ). At present, however, the conclusions drawn in 
the remainder of this work do not depend on these two 
galaxies because their metallicity errors are large. Oxy- 
gen abundances for the entire sample are listed in Table 

m 

Our inference of upper branch metallicities at M > 
2q8.2-8.5j^^ is in consistent with results reported by 
iZahid et ail (|2011h . In contrast to Figure (H Zahid et 
al. find that the mean value of R23 turns over around 
M r^ 10^-2 M0 for their z ~ 0.8 sample. They interpret 
the turnover as evidence for lower branch metallicities. 
However, as can be seen in their Figure 4, the data are in- 
creasingly noisy below M < 10®'^ M©, and objects with 
R23 > 10 have been removed as candidate AGN. The 
coupling of these effects may bias the lowest mass bin 
towards lower mean values of R23, and mimic the effect 
of a turnover in the R23 vs mass correlation. We con- 
clude that on average, galaxies at z ^ 0.6 — 0.7 fall on the 
upper branch of R23 for M > IO^^-^-SMq. At higher 
redshifts, where metallicity evolution is significant (i.e. 
lErb et al]|2006[ ). we expect that the transition from the 
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Figure 4. Our measured values of R23 decrease with increasing stellar mass, implying that most of the galaxies in our sample fall on the 
upper branch of R23. This diagram is contrasted with three other diagnostics that have previously been used to identify whether galaxies 
have low or high metallicities jMaiolino et al.ll200aiHu et al.ll200aiXia et al.ll2012l) . These include the [O III]AA4959, 5007 / H/3 flux ratio, 
the rest-frame equivalent width of H/3, and the ratio 032 = log ([O III]AA4959, 5007 / [O II] AA3726, 3729). With an H/3 EW of 352A (rest), 
object 246 from Xia et al. is not shown in the EW(H/9) panel. 

upper to lower branch of R23 will occur at higher stellar 
masses. 



5. RESULTS 
5.1. The Mass-Metallicity Relation 

Figure [5] presents the MZ relation for our 26 star- 
forming galaxies, and compares it to both the local re- 
lation (Tremonti et al. 200^, and other intermediate- 
redshift measurements from the literature. In each case, 
we have tak en care to (whe n necessary) convert stellar 
masses to a iChabrieJ (I2003D IMF (using multiplicative 
const ants given in lSavaglio et al.l [20051 iCowie fc Barged 
|2008() , and convert metall icities to a KK04 calibration 
(using equations given in iKewlev fc EllisonI 120081 ) . In 
the mass range 8.5 < Mq < 9.0, our data give a median 
value of 12 + log(0/H) = 8.63, with an RMS scatter of 
0.12 dex. 

The intermediate redshift MZ relation has been 
measured for higher mass galaxi es by sever a l au- 
thors (ISavaglio et al. 2005; Co wie fc Barged 120081 : 
iLamareille et al.l 120091 iZahid et al] 120 111) . Among the 
higher mass portion of our sample (M > 10^-^Mq, 
where we can compare directly to previous results) there 
is good agreement with all of the published MZ rela- 
tions. While the addition of our data does not distin- 



guish between these relations, we note that the Zahid 
et al. sample stands out as the largest (1350 galax- 
ies, as opposed to fewer than 100 galaxies in each of 
the o thers). These autho r s discuss the difference s be- 
tweenlSavaglio etal.l(l2005[) . ICowie fc Barged (|2008[) . and 
ILamareille et al.l (|2009l ). Thev conclude that sample se- 
lection effects and linear fits which are biased by out- 
liers may account for the differences among these works. 
While these effects are likely important, we also note 
that different d ust correction meth ods may be significant. 
On on e hand, IZahid et al.l ()2011[ ) and ILamareille et al.l 
()2009f ) calculate R23 from equivalent widths, as they are 
less sensitive to dust extinction. Although this method 
is subject to systematic uncerta inties correlated with 
galaxy colors, iLiang et"aLl (|2007| ) show that this effect 
has a relatively weak impact on metallicity (-0.2 to 0.1 
dex). Cowie & Barger, on the other hand, take the ap- 
proach of determining dust extinction from SED fits. 
However, they do not assume a higher dust extinction 
for nebular compared to stellar light, as their lower red- 
shift sample suggests that these quantities are consis- 
tent. Compa. r atively , the dust correction adopted by 
ISavaglio et al.l (|2005l ) is cruder; they use the same ex- 
tinction {Ay = 2.1) for all galaxies. If this correction 
is accurate for the median of their sample, then at high 
(low) masses the dust correction will be underestimated 
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Figure 5. The addition of our 26 low mass galaxies constrains the low mass po rtion of the MZ relation at intermediate redshifts. For 
comparison, we have re-calculated the metallicities presented in IXia et al.l II2012I ), assuming- contrary to these authors- that the three 
highest mass objects have metallicities on the upper branch. Other studies are shown for comparison, and are (when necessary) converted 
to a Chabricr (2003) initial mass function and a Kobulnickv & Kcwlcv (2004) metallicity. The dashed and dotted grey lines that follow the 
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Figure 6. The luminosity metallicity relation for the present emis- 
sion line selected sample is compared to correlations from the litera- 
ture. Notably, we compare to the measurements of the ultra-strong 
emission line galaxies reported by Hu et al. (2009). Because these 
galaxies have detected [O III] A 4363 A emission, we assume that 
they fall on the lower branch of R23. Their metallicities are recal- 
culated using the KK04 calibration (assuming no dust, as justified 
by Hu et al.). 



(overestimated), and metallicities will be overestimated 
(underestimated). Qualitatively, this effect could explain 
their steeper MZ slope. Additionally, it is worth noting 
that our data are inconsistent with an extrapolation of 



the Savaglio et al. relati on. In the sections that follow, 
we take the results from iZahid et al.l ()2011[ ) as the best 
measurement of the high mass, intermediate-redshift MZ 
relation. 

The addition of our low mass data allows for new con- 
straints on the evolution of the MZ relation. In Figure 
[5] we c ompare to the local relation measured from the 
SDSS (|Tremonti et al.ll2004D . In the mass range 10*-^ < 
M/Mq < 10^°, we find a mean metalhcity that is 0.12 
dex lower than the local relation (at Al = IO^-'^^Mq). 
This trend confirms that metallicity evolution is rela- 
tively slow from intermediate to low redshifts. Addition- 
ally, as our data are consistent wit h an extrapolation of 
the results from IZahid et al.l (|2011[ ). we agree with their 
finding that metallicity evolution is more significant at 
lower masses than at higher masses. This trend is qual- 
itatively consistent with downsizing in the later phases 
(z ~ 1 — 2 to today) of galaxy evolution. After high mass 
galaxies have assembled most of their stars they have also 
made most of the metals that they will ever produce. At 
the same time, these higher mass galaxies are less ef- 
fective at removing metals throu gh winds, and they are 
also less efficient at accreting gas (jDekel et al.ll2009[ ). As 
a result, higher mass galaxies are less capable of changing 
their metallicities after they've assembled most of their 
stars. Under this scenario, we would expect lower mass 
galaxies to evolve onto the present MZ relation at later 
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Table 2 

Emission Line Measurements 



IMACS ID 


COSMOS ID 


RA 
(J2000) 


Dec 
(J2000) 


2, 


F([0 III] A5007) 


F([0 III] A4959) 


F(H/3) 


F([0 II] A3727) 


EW{H/3) 
A (rest) 


27.5+5-0.47 


799604 


10:00:11.314 


+02:04:07.60 


0.675 


16.7 ±6.3 




39.2 ±0.7 


53.9 ±1.3 


12.7 


99.5-3-0.84 


1235118 


10:00:53.842 


+02:22:06.88 


0.685 


19.75pml2.7 




25.1 ±0.6 


58.0 ±1.2 


8.5 


90.5-1-6-0.63 


1265343 


10:00:09.101 


+02:19:52.02 


0.678 


10.2 ±2.7 




6.3 ±0.3 


16.4 ±0.4 


9.2 


100.5-2-0.41 


1234271 


10:00:50.551 


+02:22:22.66 


0.683 


13.4 ±2.2 




19.0 ±0.4 


48.5 ±0.9 


9.7 


65.5-1-0.69 


1003946 


10:00:47.710 


+02:13:36.83 


0.620 


26.0 ±0.4 




19.4 ±3.0 


40.9 ±0.7 


9.6 


87.5-1-4-0.86 


1242162 


10:00:18.502 


+02:19:06.82 


0.656 


15.4 ±0.7 


5.4 ±0.4 


5.4 ±1.2 


19.4 ±0.8 


9.2 


95.5-1-4-0.94 


1236769 


10:00:21.331 


+02:21:05.22 


0.667 


158.5 ±1.4 


52.5 ±0.5 


40.8 ±0.5 


112.3 ±1.0 


32.0 


35.5-1-5-0.23 


794540 


10:00:09.038 


+02:06:08.16 


0.620 


57.0 ±0.4 


19.1 ±0.4 


19.6 ±0.5 


70.5 ±1.0 


14.8 




1015084 


10:00:41.863 


+02:09:12.66 


0.678 


309.7 ±0.9 


104.1 ±0.5 


75.1 ±0.4 


200.3 ±1.0 


39.5 


54.5-1-6-0.09 


1036285 


10:00:01.553 


+02:10:51.89 


0.684 


53.3 ±1.0 


17.2 ±0.6 


18.5 ±0.3 


49.6 ±0.6 


27.0 


25.5-1-7-0.08 


801017 


09:59:58.039 


+02:03:38.89 


0.638 


36.9 ±0.4 


12.8 ±0.4 


13.2 ±0.8 


46.9 ±1.0 


12.1 


23.5+2-0.33 


775548 


10:00:32.287 


+02:03:09.25 


0.670 


32.8 ±0.7 


12.2 ±1.3 


14.8 ±0.5 


44.8 ±1.1 


22.3 




1035959 


09:59:59.347 


+02:10:57.10 


0.628 


35.2 ±0.3 


12.7 ±0.2 


13.2 ±0.7 


34.7 ±0.6 


11.7 




982000 


10:01:07.140 


+02:12:08.99 


0.621 


235.1 ±0.6 


77.2 ±0.5 


53.3 ±0.7 


112.7 ±0.9 


44.0 




1013478 


10:00:18.701 


+02:09:41.74 


0.635 


34.9 ±0.3 


12.1 ±0.3 


11.2 ±0.5 


40.7 ±0.6 


28.4 




980802 


10:01:07.135 


+02:12:38.22 


0.621 


27.3 ±0.4 


9.2 ±0.6 


11.9 ±0.7 


40.5 ±0.7 


49.3 


28.0+7.1=^ 


799190 


10:00:03.035 


+02:04:13.98 


0.638 


167.5 ±0.6 


57.3 ±0.4 


35.2 ±0.8 


91.3 ±1.0 


56.3 




770439 


10:00:42.334 


+02:05:43.53 


0.629 


34.4 ±0.4 


11.5 ±0.4 


9.8 ±0.7 


23.9 ±0.8 


22.9 




1009808 


10:00:44.558 


+02:11:18.04 


0.616 


96.9 ±0.4 


32.5 ±0.5 


22.2 ±0.9 


45.2 ±0.8 


38.0 




1237667 


10:00:52.368 


+02:20:58.18 


0.640 


59.2 ±0.4 


21.2 ±0.4 


13.6 ±0.4 


24.9 ±1.0 


40.9 




772773 


10:00:25.567 


+02:04:44.36 


0.616 


21.5 ±0.5 


8.2 ±0.8 


6.9 ±1.0 


22.3 ±1.0 


11.3 


48.5-4-0.66 


989316 


10:01:08.537 


+02:09:21.70 


0.674 


13.8 ±0.5 


4.7 ±1.0 


4.96 ±0.2 


16.1 ±0.7 


17.7 


54.5+4-0.18 


1010994 


10:00:16.097 


+02:10:51.60 


0.634 


28.7 ±0.2 


10.0 ±0.3 


5.4 ±0.5 


13.4 ±0.5 


12.2 


47.5+3-0.85 


1015516 


10:00:27.792 


+02:09:06.43 


0.679 


34.1 ±0.3 


12.4 ±0.2 


8.5 ±0.2 


17.7 ±0.6 


43.5 




1010583 


10:00:21.408 


+02:11:00.43 


0.639 


45.2 ±1.5 


14.6 ±0.8 


10.6 ±2.1 


28.5 ±0.6 


27.9 




1217842 


10:00:54.658 


+02:19:01.81 


0.633 


21.4 ±0.3 


8.4 ±0.4 


5.1 ±1.7 


9.8 ±0.8 


19.6 






Objects lacking clear optical counterparts in 


the COSMOS data 








52.5+8-0.69 




09:59:51.374 


+02:10:23.28 


0.678 


15.6 ±0.8 


3.7 ±0.9 


2.1 ±0.1 


3.6 ±0.4 




45.5+7-0.25 




09:59:55.699 


+02:08:38.51 


0.678 




2.2 ±0.7 


2.6 ±0.4 


9.6 ±1.0 


6.5 


46.5-7-0.04 




10:01:25.152 


+02:08:52.34 


0.672 


6.2 ±0.5 


2.5 ±1.5 


2.4 ±0.8 


5.4 ±0.7 




65.5-2-0.14 




10:00:50.904 


+02:13:36.52 


0.624 


14.1 ±0.3 


5.0 ±0.9 


5.4 ±0.2 


6.8 ±0.8 




47.5+6-0.90 




10:00:07.145 


+02:09:07.81 


0.639 


14.4 ±0.4 


5.0 ±3.8 


3.4 ±0.8 


<9.5 




87.5-1-0.06 




10:00:48.792 


+02:19:06.54 


0.620 


15.4 ±0.3 


5.3 ±1.6 


6.0 ±1.7 


1.9 ±1.3 




Candidate AGN 


55.5-3-0.60 


984577 


10:01:01.126 


+02:11:07.85 


0.623 


39.5 ±0.4 


13.7 ±0.4 


16.4 ±0.5 


51.6 ±0.7 


9.9 


73.5+3-0.54 


998509 


10:00:25.046 


+02:15:38.15 


0.642 


38.7 ±0.5 


12.4 ±0.6 


48.4 ±0.8 


93.3 ±0.9 


9.4 



Note. — Emission line fluxes are in units of 10 erg s ^ cm . Equivalent widths and H/9 fluxes are given as measured (without applying 
the lA stellar absorption correction.) Objects lacking an IMACS ID were selected strictly from the COSMOS narrowband (NB816) imaging, as we 
described in ^ The [O II] flux is the sum of the AA3726, 3729 doublet fluxes. In six cases, one of the lines in the [O III] doublet lines fell on an OH 
sky line; in this situation, metallicities are derived by assuming a flux ratio of F([0 III] A5007)/ F([0 III] A4959) = 3. 

'^ This galaxy was drawn from our earlier (wider and shallower) survey first presented in Marti rTet al.l 120081) . 



cosmic times. Hence, extending observations to lower 
masses at all redshilts will provide essential insights into 
the physics of galaxy assembly and metal production. 
We compare our data to theoretical predictions for the 
IVIZ relation in gT] 

5.2. The Luminosity- Metallicity Relation 

In Figure ini we show the metallicity luminosity relation 
derived from our data. While the IMZ relation may repre- 
sent a more fundamental correlation than the metallicity- 
luminosity relation, the latter allows comparison to stud- 
ies where stellar-mass measurements have not been pos- 
sible. Here, we focus on the ultra-strong e mission line 
(USEL) sample presented by IKakazu et al.l f2007) and 
iHu et all ()2009D . Like the galaxies in our sample, the 
USELs lie at intermediate redshifts and are emission line 
selected. However, the majority of the USELs are low 
metallicity galaxies with detected [O III] 4363 A emis- 
sion. In Figure [HI we compare the metallicity-luminosity 
relation of the USELs to our sample. In order to make 
a fair comparison, we have calculated the metallicities 



of the USELs using R23 and KK04, assuming lower- 
branch so lutions. No correction for dust is applied, which 
IHu et al.l p009) argue is appropriate for these faint, low 
mass objects. It is worth noting that, under these as- 
sumptions, the majority of USELS have 12 -|- log (0/H) 
> 8.0, and would not be classified as ex tremely metal 
poor galaxies (12 + log(0/H) < 7.65 ; iKniazev et alj 
|2003( ). However, this classification changes when direct 
method metallicities are used: in this c ase, the oxygen 
abundances reported bv lHu et al] (|2OO90 are 0.03 - 0.74 
dex lower (0.4 dex on average). 

Figure H] shows clearly that the USELs have lower 
metallicities than the emission-line selected galaxies in 
the present work. In fact, this trend is not surpris- 
ing. Because of the differing selections for the USELs 
and the present sample, the former have H/? equiva- 
lent widths that are two times larger (on average) than 
the latter. This difference translates to higher specific 
(mass-normalized) SFRs for the USELs. One possible 
explanation is that more actively star-forming galaxies 
will have B-band luminosities which are high for their 



Table 3 

Derived Properties 
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IMACS ID COSMOS ID Mg U - B 



log M* 



log SFR 



iV 



log R23 



log 032 12+log(0/H) 



log q 



27.5+5-0.47 
99.5-3-0.84 
90.5+6-0.63 
100.5-2-0.41 
65.5-1-0.69 
87.5+4-0.86 
95.5+4-0.94 
35.5+5-0.23 

54.5+6-0.09 
25.5+7-0.08 
23.5+2-0.33 



28.0+7.1 



48.5-4-0.66 
54.5+4-0.18 
47.5+3-0.85 



799604 


-20.93 


0.52 


10.0510:00 


1235118 


-20.71 


0.45 


9-77lH^ 


1265343 


-19.24 


0.61 


q 74+0.00 


1234271 


-20.39 


0.42 


9 62+0-08 


1003946 


-19.56 


0.42 


q 17+0.09 


1242162 


-19.15 


0.44 


q 1K-O.O6 
y-J^O-0.16 


1236769 


-20.08 


0.37 


9-121H? 


794540 


-19.89 


0.30 


9.O2IHI 


1015084 


-20.26 


0.21 


8.99lHg 


1036285 


-18.93 


0.40 


8.98j:0;i0 


801017 


-19.31 


0.37 


Q7+0.05 


775548 


-19.32 


0.34 


« Sfi+0-06 


1035959 


-19.11 


0.37 


oq+0.09 


982000 


-19.54 


0.25 


og + O.OS 


1013478 


-18.64 


0.36 


79+0.08 


980802 


-18.74 


0.37 


cq+0.09 


799190 


-19.07 


0.24 


cn+O.OO 
8.62_o.io 


770439 


-18.46 


0.36 


c-q+O.lO 


1009808 


-18.62 


0.28 


8 45+0-18 


1237667 


-17.99 


0.30 


^q+O.Ol 


772773 


-18.72 


0.28 


8.40t0-J0 


989316 


-17.94 


0.35 


8.39t»;f3 


1010994 


-18.18 


0.40 


,0+0.07 


1015516 


-17.82 


0.36 


,7+0.73 


1010583 


-18.00 


0.26 


99+0.20 
»-^^-0.16 


1217842 


-18.00 


0.20 


7 Qo+0.09 
7.98 „ (If; 



1.23 

0.95 

-0.54 



+0.34 
-0.21 
+0.20 



+0.12 
-0.00 
n ^9+0.38 

r + 0.44 

-0.06 
+0.07 

0.00 
+0.59 

0.22 



-0.45 

-0.42 

0.30 



09+0-22 
'-'-'^^-0.05 
Q 33+0.15 
'-'-'-'•-'-0.00 



+0.38 

0.46 
+0.00 

0.25 
+0.18 

0.16 
+0.20 

0.18 
+0.23 

0.48 
+0.54 

0.15 
+0.33 

0.50 
+0.46 

0.20 
+0.46 

0.00 
+0.18 

0.35 
+0.76 

0.02 
+0.30 

0.02 
+0.37 

0.52 
+0.06 

0.01 

j+0.17 
0.32 



0.08 
-0.06 
-0.33 
-0.36 

0.43 
-0.28 
-0.56 
-0.04 
-1.08 
-0.13 
-0.51 
-0.61 
-0.21 
-1.41 



-1+0.13 
-0.48 

-1+0.08 
-0.19 

-1+0.30 
-0.05 

1+0.28 
-0.31 

-1+0.07 
-0.21 
f0.18 
-0.00 

1+0.27 
-0.40 

20+0-40 
'^•^'^-0.07 

00+0-00 
'^•'-"J-0.00 

+0.44 

0.60 

+0.22 



I.8O: 
1.40J 
0.20J 
l.OOJ 
0.60J 

o.ooj 

0.60^ 



0.60 
0.20 
0.20 
0.20 
0.40 
0.60 
0.20 
0.20 
0.20 
0.20 



0.20 

+0.14 
0.20 

+0.23 
0.20 

+0.24 
0.26 
0.46 
0.24 

+0.27 
0.20 
0.43 
0.07 

+0.11 
0.20 

+0.27 
0.20 



n nn+o-82 

U.UU_g QQ 
QQQ + 0.15 

'J-^'J-0.00 



-0.48] 



0.60 
0.00 
1.60 
0.60 
0.00 



+0.40 
0.60 
0.04 
0.00 

+0.30 
1.54 

+0.38 
-0.56 

+0.22 



0.68 + 0.10 
0.77 + 0.07 
0.63 + 0.08 
0.71 + 0.10 
0.62 + 0.09 
0.78 + 0.12 
0.85 + 0.05 
0.85 + 0.06 
0.89 + 0.01 
0.87 + 0.08 
0.83 + 0.06 
0.78 + 0.04 
0.75 + 0.06 
0.92 + 0.03 
0.96 + 0.08 
0.83 + 0.05 
0.96 + 0.03 
0.84 + 0.04 
0.89 + 0.03 
0.87 + 0.04 
0.81 + 0.08 
0.89 + 0.09 
0.92 + 0.06 
1.06 + 0.12 
0.96 + 0.11 
0.85 + 0.16 



-0.94 + 0.16 

-0.78 + 0.22 

-0.14 + 0.11 

-0.75 + 0.19 

-0.26 + 0.05 

0.03 + 0.04 

0.09 + 0.10 

-0.03 + 0.07 

0.32 + 0.02 

-0.03 + 0.14 

-0.04 + 0.06 

-0.06 + 0.05 

0.08 + 0.06 

0.32 + 0.08 

-0.12 + 0.11 

-0.11 + 0.07 

0.33 + 0.07 

0.22 + 0.04 

0.39 + 0.07 

0.51 + 0.08 

0.12 + 0.04 

-0.13 + 0.14 

0.46 + 0.02 

-0.08 + 0.25 

0.13 + 0.13 

0.49 + 0.05 



1+0.08 
-0.12 

O fiq+0.08 

».D3_o 10 

a+0.07 
'-0.06 
,+0.13 

-0.10 
1+0.09 

-0.06 
j+0.18 

-0.10 
,+0.16 

-0.10 
j+0.09 

-0.06 
7+0.01 

-0.02 
7+0.02 

-0.13 
1+0.10 

-0.08 
,+0.05 

-0.04 
;+0.07 

-0.06 
,+0.08 

-0.04 
;+0.21 

-0.16 
;+0.10 

-0.07 
j+0.07 

-0.10 
1+0.07 

-0.08 
7+0.07 

-0.04 
,+0.07 

-0.06 
■,-f0.15 
" -0.06 
a+0.02 

-0.18 
1+0.08 

-0.09 



;.80"* 



8.i 

8.76:' 

8.89:1 

8.68^ 

8.41^ 

8.58:* 

8.57J 

8.57J 

8.61^ 

8.72^ 

8.75:! 

8.51:1 

8.35^ 

8.65^ 

8.43:1 

8.64J 

8.57I 

8.63:! 

8.66:! 

8.49^ 

8.49"* 



^-40.0. 30 

-^«-0.26 

o An+0.50 
8.4U_o,40 



7.05] 

7.65 



7 00+0-00 
'•'^^-0.13 

-+0.14 

0.10 

+0.12 

-0.07 

7 -I -I +0.06 

'-11-0.10 

7 c:c:+0.09 

'■''•^-0.11 

7 66+0-" 
'■°°-0.12 



7.57 



+0.14 
-0.12 



7 56+0-07 
'-^^"-0.12 

7 oq+0.02 
'■"■-'-0.03 
+0.14 
0.05 
+0.09 
0.08 
+0.09 
0.08 



7.54 

7.56 
7.61 

7.74 



+0.10 
-0.11 



7.80«:lg 

7A0t°f, 
+0.11 

0.07 
+0.07 
0.13 
+0.06 
-0.07 
+0.11 
-0.07 

8.05_(,:27 

7 79+0.08 
'■'^-0.06 
+0.03 
0.18 
+0.05 
0.05 
+0.22 
-0.19 
+0.15 



7.53 
7.77 
7.80 
7.91 



7.44 
7.91 
7.36 
7.60 
8.05 



+0.22 



Objeets with ambiguous eounterparts in the COSMOS data 



52.5+8-0.69 
45.5+7-0.25 
46.5-7-0.04 
65.5-2-0.14 
47.5+6-0.90 
87.5-1-0.06 




1.03 + 0.04 
0.74 + 0.12 
0.77 + 0.16 
0.68 + 0.03 
0.86 + 0.14 
0.57 + 0.13 



0.73 + 0.05 
-0.04 + 0.13 
0.21 + 0.10 
0.45 + 0.06 
0.61 + 0.35 
1.04 + 0.29 



Candidate AGN 



55.5-3-0.60 
73.5+3-0.54 



984577 -19.78 
998509 -21.23 



0.62 
0.54 



10.25 



+0.01 
0.00 



10 23+000 



0.67: 
0.62 



1.20 



+0.00 

0.16 
+0.55 -, W 

0.08 ^^^"-0.12 



-0.10 
-0.20 
-0.20 



Note. — These quantities are derived from COSMOS imaging and our measured line ratios. Magnitudes and colors are in AB units and are 
derived by integrating the best-fitting SED under the appropriate bandpasses. Stellar masses and SFRs are giv en in solar units, and were derived 
using FAST iKriek^et al. 2009). For this derivation we used a Chabrier (2003) initial mass function and a Calz etti et al.l 1I2OOO1 ) extinction curve. 
The line ratios, metallicity, and ionization parameter (q) are derived, taking into account a lA (equivalent width) correcti on to the H/3 fluxes, as 
well a s the dust correction and its associated uncertainty. Metallicities and ionization parameters are calculated using the IKobulnickv &: Kewlevl 
H2004I ) calibration. A dust correction has not been applied to the objects have ambiguous optical counterparts in the COSMOS imaging survey. 
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stellar mass. In this scenario, galaxies that fall along 
the MZ relation will move to the right when plotted 
in the luminosity-metallicity relation, as in Figure [6l 
However, another explanation is possible. Within the 
local MZ relation, galaxies with higher specific SFRs 
have lower metallicities and follow a mass-metallicity- 
SFR plane c alled the Fundamental Metallicity Re- 
lation (FMR: | Lara-Lopez et al.l 120101: [M annucci et al.l 
[2010l[20lH [Yates et al.l(2012l:IAndre\TOfc Martini 201^7 
Hence, the differing metallicities between the USELs and 
the present sample supports the existence of the FMR in 
the relatively unexplored low mass, intermediate redshift 
regime. We next provide the first quantitative test of the 
FMR for intermediate redshift galaxies below 1O^M0. 



6. THE FUNDAMENTAL METALLICITY RELATION 

In the previous section, we showed that, for a fixed lu- 
minosity (and possibly mass), galaxies with higher equiv- 
alent width emission lines (i. e. the USELs) have lower 
metallicities. This result supports the idea that star for- 
mation drives some of the scatter in the MZ relation, 
even out to intermediate redshifts. For local galaxies 
at fixed stellar mass, objects with highe r SFRs have 
lower metallicities than lower- S FR ga l axies (lEUison et al.l 
2008; Mannucci et al. 20l dl |201H lYates et al.l 120121: 
Andrews fc Martinil 120131) . iMannucci et al.l (|2010D de- 
scribe this relation as a plane, and refer to it as the 



Fundamental Metallicity Relation (FMR). Remarkably, 
Mannucci et al. report that high-redshift galaxies also 
lie on the plane (albeit at lower masses and higher SFRs 
than most low-z galaxies), implying that the physics 
governing the evolution of metal enrichment has not 
changed over cosmic time. This finding has recently 
been confirmed for a large sample of intermedia te red- 
shift (z ~ 0.6) galaxies from zCOSMOS (Cres ci et al.l 
l2012f ). Nevertheless, tests for evolution of the FMR have 
not yet been extended below M '^ 10^ M©. Here, we 
provide the first such test. 

In light of the significant differences between metallici- 
ties d erived from various calibrations (jKewlev fc EUisoiil 
l2008f ). it is essential that we use metallicities which are 
derived consistently with those in (Mannucci et al. 2010, 
2011; the latter work extends the initial analysis to lower 
masses). These metallicities a re based on the calibra- 
tion from lMaiolino et al.l ()2008[ ). which is semi-empirical. 
The high metalli cities are derived from the photoion- 
ization models of iKewlev fc Dopi"tal ()2002t ). whereas the 
lower metallicities are constrained from direct T^ mea- 
surements reported by Nagao ct al. (2006). While the 
Maiolino et al. calibration of R23 does not include a de- 
pendence on the ionization parameter, in other calibra- 
tions (i.e. KK04) this dependence is weak on the upper 
branch of R23, where the present sample lies. There- 
fore, for this section, we re-calculate the met allicities for 
our sample using the iMaiolino et al.l (|2008[ ) calibration 
of R23. 

Figure [71 compares our galaxies to the z = 0.1 FMR re- 
ported in iMannucci et al] (|201lD . On the top panel, we 
show the MZ relation, with our data points color coded 
by three different SFR bins. The FMR, which is evalu- 
ated and shown for the same three SFR bins, is given by 




1.41 < log SFR < -0.42 

0.36 < log SFR < 0.02 

— 0.08 < log SFR < 1.23 



9.0 9.5 

Log M/M „ 



Figure 7. The fund amental metallicity relation, derived using the 
IMaiolino et al.l 1)200 8') metallicity calibration, shows that low mass 
galaxies with high SFRs have systematically lower metallicities. 
The galaxy 23.5-1-2-0.33, which strongly influences the weighted 
mean of the residuals (see text), is shown as an open green sym- 
bol. Top- The MZ relation is plotted, with galaxies color-coded 
in three bin s of SFR. T he curves show the local FMR reported by 
[Mannucci et al.l II2011I '). FMR. Bottom- Residuals from the local 
FMR (data minus relation given in Equation [TJ show increased 
scatter towards low stellar masses, owing to larger measurement 
errors for faint galaxies. (Note that residuals are calculated using 
measured SFRs rather than the mean values where Equation [T] has 
been evaluated for the top panel. As a result the residuals in the 
bottom panel do not exactly correspond to the differences between 
the data and models in the top panel.) Inset- Likelihood functions 
for the mean of the residuals are shown for different subsets of the 
data, including the full sample, low and high mass subsamples, and 
the SFR divided subsamples. With one exception, the likelihood 
functions shown here include the galaxy 23.5+2-0.33, which drives 
some of the trends towards larger positive residuals (see text). The 
thin green curve shows the effect of removing this outlier. 



the relation (Equation 2 from IMannucci et al.ll201ll) : 

12 + log(0/H) = 8.90 -h 0.37m - 0.14s - 0.19™^ 

+0.12ms - 0.054s2 for /10.32 > 9.5 

= 8.93 -I- 0.51(^0.32 - 10) for /io.32 < 9.5, 

(1) 
where m = log(Af) — 10, s = log(SFR), and /i32 = 
log(M) — 0.321og(SFR). This exercise shows clearly that 
a planar relation exists to lower masses at 2; '^0.6 — 0.7, 
with low SFR galaxies showing higher metallicities. To 
highlight the differences from the FMR, the bottom panel 
shows the residuals from the local plane. Here, we have 
included the uncertainty in the SED-derived SFR, as it 
enters into the metallicity that is predicted by the local 
FMR. Overall, we see that most galaxies have metallici- 
ties that fall within 0.2 dex of the local FMR, However, 
at low masses and low SFRs there are more objects with 
Alog (0/H) > 0, so it is reasonable to question whether 
these galaxies deviate slightly from the FMR. There- 
fore, we calculate the mean of the residuals, /i = (Alog 
(0/H)), and its uncertainty by using a maximum likeli- 
hood estimation to account for the individual measure- 
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Table 4 

Residuals from the local FMR 


Sample 




Mean of the Residuals, /x 


Full sample 
M > lO'^-S Mq 
M < lO^-S Mq 
-1.41 < log SFR < 
-0.36 < log SFR < 
0.08 < log SFR < 


-0.42 
0.02 
1.23 




0.016 + 0.020 
-0.074 ±0.083 
0.046 ± 0.022 
0.013 ±0.035 
0.086 ± 0.033 
-0.042 ±0.029 


Without 23.5+2-0.33 


Full Sample 

-0.36 < log SFR < 0.02 

M < lO^-S Mq 




-0.009 ±0.021 
0.027 ±0.043 
0.018 ±0.025 



ment errors. Under this approach, the probability of 
observing a galaxy with a residual yi and error, ai is 
Pi = e"(^-"^)'/2<Ti^^ rpj^g^^ jpj. ^j^g j^jj sample the likeli- 
hood function becomes C{fj,) — Y[Pi- We evaluate this 
likelihood function between —0.4 < /i < 0.4, taking the 
maximum and the 68% confidence intervals of C{fi) to 
represent the mean and its error. This estimate is made 
for the full sample, as well as low and high mass subsam- 
ples (divided at 10^-^ Mq) and the subsamples in three 
bins of SFR. 

In the inset panel of Figure [71 we show the likelihood 
functions for the mean of the FMR residuals. These cal- 
culations indicate only insignificant deviations from the 
local FMR. In Table |4j we give the mean of the residuals 
and its uncertainty for the full sample and various sub- 
samples. At first glance, we see two weakly significant 
deviations from the local plane: one among the lower 
mass {M < 10^'^ M0) galaxies, and the other for those 
with intermediate SFRs (green points). However, upon 
closer inspection we find that these trends arise solely 
because of one galaxy (23.5+2-0.33) with a metallicity 
3.4cr above the fundamental plane. (Furthermore, this 
galaxy has a strong influence on our calculation not be- 
cause its residual is large, but rather because its errors 
are small.) Without this galaxy, the deviation disappears 
(see the thin green curve in the inset of Figure [7]) , and 
the mean of the residuals fall within l.Scr of zero for all 
of the subsamples listed in Table U) Since any claims 
of a deviation from the FMR should be based on more 
than one object, we conclude that our data show no com- 
pelhng evidence for evolution at 10^ < M/Mq < 10^-^. 
This result lends greater leverage to the finding that the 
FMR is not evolving for hi gher mass galaxies at these 
redshifts ()Cresci et al.ll2012D . 

It is also interesting to examine the scatter about 
the FMR, as intrinsic scatter could indicate that there 
are additional physical processes that infiuence gas- 
phase metal li cities . For intermediate redshift galaxies, 
iCresci et al.l (|2012f ) report a scatter of 0.11 to 0.14 dex 
about the FMR. Yet, in our sample, the scatter appears 
to be larger at lower masses. We find an RMS of 0.2 dex 
below M < 10^'^ M0. Assigning some galaxies to the 
lower branch of R23 does not reduce this scatter. The 
amplitude of the residuals are larger (or the about the 
same for a few galaxies) when the lower metallicity solu- 
tions are adopted. In order to disentangle intrinsic scat- 
ter from increased measurement errors in faint galaxies, 
we use a Monte Carlo simulation. First, we assume that 



there is no intrinsic scatter relative to the FMR. Then, 
focusing only on the galaxies with M < 10^'^ M©, we 
generate 10,000 mock realizations of the FMR residuals, 
perturbing the data in Figure 7 by their uncertainties. 
We calculate the RMS about the FMR for each real- 
ization, and find a mean of 0.23 dex, in good agreement 
with observations. This exercise implies that the amount 
of scatter that we observe can be entirely explained by 
our measurement errors. Therefore, we conclude that 
the intrinsic scatter about the FMR must be small. Re- 
peating the simulation with additional intrinsic scatter 
shows allows us to place an upper limit on this scatter. 
We find that an intrinsic scatter of 0.16 dex produces an 
observed RMS < 0.2 dex in only 5% of Monte Carlo re- 
alizations. Remarkably, the scatter about the FMR for 
the present sample is significantly smaller than the 0.4 
dex RMS found for similar mass galaxies in the local uni- 
verse f|Mannucci et al. 20111) . We therefore conclude that 
the increased scatter in the low mass, z = 0.1 FMR is 
likely the result of increased measurement uncertainties 
for faint galaxies. 

7. CONSTRAINING GAS FLOWS WITH THE 
MASS-METALLICITY AND SFR-M* RELATIONS 

As we introduced in §1, the MZ relation is an im- 
portant measure of galaxy evolution. Galaxies are not 
closed boxes; gas must be accreted from the IGM t o 
fuel star formation (e.g. Erb 2008; Bouche et al.l 120101) , 
and galactic outfiows are commonl y observed at both 
low and high- redshift (e . g. Heckman et al.l Il990t IMartin 
20051: iHenrv et al.l 120071: iSoto et al.l 120121: iKornei et al. 



2012[ lErb et al.ll2012D . The MZ relation provides an in- 



dependent, albeit indirect, constraint on these gaseous 
flows. On one hand, the accretion of primordial gas will 
reduce the gas-phase metallicity in a galaxy. On the 
other hand, galactic winds can create a deviation from 
the closed-box model by low ering the gas fraction with- 
out c hanging the metallicity (|Edmundall990l : iDalcantoiil 
l2007f) . Combining our data with the higher mass MZ re- 
lation provides an important increase in dynamic range 
that allows for more stringent constraints on models. 

Our emission line selection of galaxies does not strongly 
bias our sample or hamper our ability to test theoretical 
models. In Figure [S] we show that our sample is consis- 
tent with the DEEP2 MZ relation reported by Zahid ct 
al. (2011; at least in the limited mass range where our 
data overlap). This agreement can be understood, since 
the FMR dependence on the SFR is relatively weak. For 
example, if the SFRs of our galaxies wer e biased high by 
a fac tor of several, then (following the iMannucci et al.l 
I2OIII FMR) an "unbiased" sample would have metallici- 
ties which are approximately 0.1 dex higher. An increase 
in metallicity this large is implausible, as it would imply a 
lack of metallicity evolution inconsistent with other stud- 
ies. Likewise, it is worth pointing out that the high S/N 
of the spectra shown in Figures [T] and [5] implies that our 
sample is not biased towards higher metallicities by a 
requirement for strong H/3 emission. The objects which 
were excluded due to low S/N all fell among the class 
#3 objects discussed in §2.1; these are either extremely 
faint, low- mass galaxies with uncertain optical counter- 
parts, or brighter galaxies that fell partially outside our 
slit. These galaxies have poorer quality spectra overall; 
it is unlikely that their exclusion biases our results. 
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Figure 8. The SFR-M, correlation (top) and MZ relation (bottom) constrain models from lDave et all 1)20121 ). Large red squares show 



the average mass, metallicity, and SFR for four mass bins: M > 10^-^Mq, 10**--' < M/Mq < 10^ 



10» 



< M/Mq < 10«-y and M < 10« 



M0. Error bars on the binned data represent the standard deviation of the mean. We combine our data with binned higher mass data 
from DEEP2: the z = 0.5 - 0.7 SFR - M, correlation reported by Noeske et al. (2007; converted to a Chabricr 2003 IMF), and the z ~ 0.8 
MZ relation from lZahid et aP 1)20111 ). MZ models are normalized by taking a nucleosynthetic yield that minimiz es v ^ for the union of our 
data and the Zahid et al. data. SFR-M* models are no t renormalized. The models are described in detail in i|7.ll The halo mass scale 
shown at the top assumes the M« — Mf^^lo relation from lMoster et aLl 1)20101) . 
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In order to further verify that our sample is not biased 
by emission-hne selection, we place our galaxies on the 
SFR-M* relation (i.e. the star-forming main sequence) 
in the top panel of Figure [51 This comparison demon- 
strates that our galaxies do not have extreme SFRs for 
their stellar masses. Rather, they lie on an extrap olation 
of the relation reported by iNoeske et al.l ()2007() . Our 
sample could only be significantly biased if the "true" 
SFR-M* relation turns over and becomes steeper at low 
stellar masses. While this trend is predicted to occur 
when photoionization heating suppresses star formation 
in low mass halos, it is not expected to be im portant 
above Mhaio > lO^M© (jOkamoto et all [20081) . While 
beyond the scope of this paper, it is worth noting that 
we do not detect the flattening of the SFR-M^, r elatio n 
at low masses which is reported by iPirzkal et al.l (|2012l ) . 

7.1. Theoretical Framework 

One set of models which have recently gained traction 
are those that adopt an equilibrium between mass inflow 
(Min), outflow (Mout), and gas consumption via star 
formation (M.: iBouche et al.l 120101: iDutton et al.l 120101: 



iFinlator fc Dav^l2008t iDave et al.ll2012D . These models 
provide a relatively simple interpretation of the FMR as 
a manifestation of stochastic star formation. When a 
galaxy accretes gas, its metallicity is diluted and its star 
formation increases. The increased star formation serves 
to consume the extra gas and return the galaxy to its 
equilibrium. Likewise, a pause in gas accretion will lower 
star formation rates and increase gas-phase metallicities 
until accretion resumes. These models not only predict a 
MZ relation, but also the SFR-M* relation that we have 
shown in Figure [8l above. Therefore, in the remainder 
of this section we leverage the two relations to gain the 
best possible constraints on the models. 

Quantitatively, we can express the eq uilibr ium cond i- 
tion as (reproducing Equation 1 from lDave et al.ll2012[ ): 

mn = Mout + M^- (2) 

Following this ass umption, IFinlator fc Pavel (|2008| ) and 
IDave et akl (|2012[ ) show that SFRs and metallicities can 
be written as: 



SFR^ 



CM. 



grav 



and 



Zn = 



(l + 7?)(l-a,) 

_y i__ 

1 + r]l — az 



(3) 



(4) 



The quantities in Equations [3] and H are defined as fol- 
lows: 

1. Zq is taken to represent the mass fraction of oxygen 
in the ISM. This quantity is converted to the units 
of 12 + log(0/H) by taking log(Zo) = log(0/H) 
-I- log( I X Mo/Mh ). Mo and Mh are the atomic 
masses of oxygen and hydrogen. 

2. Mgj-av is t he cosmologi c al bar yonic accretion rate, 
taken from lDekel et al.l (|2009D . Because this quan- 
tity (as well as others below) are expressed in terms 
of halo mass, we ado pt the M^—Mhaio relation from 
iMoster et al.l(|20lol ). 



3. C is a quantity that represents preventive (rather 
than ejective) feedback from gas heating. There 
are multiple forms of preventive feedback, and they 
combine multiplicatively. In the range of halo- 
masses where our objects lie, the most important 
of these is gravitational heating from virial shock 
forma tion in accreted gas. Following IDave et"al] 
(l2012f). we take the analy tical form given by 
iFaucher-Giguere et all (|201lD : 



Q 



0.47 



M, 



halo 



1012 A/, 



(5) 



At high masses star formation quenching associ- 
ated with supermassive black holes becomes impor- 
tant. While we do not ai m to constrain quen ching, 
we implement it following IDave et al.l (|2012l ): 



^quench 




-1.5 



(6) 



The mass scale that we have adopted for Cquench 
is for z = 0, and while it may b e higher at earlier 
cosmic times (iDekel et al.l l2009l. the exact scaling 
is unimportant in the present analysis. At the low 
mass end, we take the photoionization feedback to 
be unimportant; as mentioned above, it is not ex- 
pected to play a significan t role at Mhaio ^ 10^ Mq. 
Finally, IDave et"aLl (|2012D do note that heating by 
galactic winds may be important, but the physics 
of this effect is poorly understood. As such, they 
adopt an arbitrary parameterization of Cwinds to 
consider its qualitative effects. As we will show 
below, this additional heating is not needed to re- 
produce the present data. 

4. y is the nucleosynthetic yield o f oxygen, which 
is between 0.008 < y < 0.021 (jFinlator k Davel 
[20081. Since metallicity calibrations are uncertain 
(jKewlev fc Ellisonll2008D . we take y as a free pa- 
rameter as we search for a normalization that best 
matches our observed MZ relation. Hence, our 
comparison to the MZ relation concerns its slope, 
but not the normalization. 



77 is the mass loading factor, which is defined 
as the^ ratio of the outflow rate to the SFR: 
77 = Mout/M^,. In the comparison that follows, 
we adopt mass-loading factors that correspond 
(respectively) to momentum- and ener gy- driven 
wind s (i.e. [M urray ct al. 2005; Finla tor fc Davel 
20081: IDave et all 1201 Ibt IPeeoles fc Shankari 120111: 



Dave et al.ll2'oTa : 



V = (Mhalo/IO^^MQ)-^/^ 

11 = {Mhaio /iQ^^MQy^'\ 



(7) 
(8) 



The norma l izatio ns of these 77 are taken from 
IDave et al] (|2012D . In addition to these two pa- 
rameterizations of 77, we compare to the case of no 
winds (77 = 0), also 77 = 1, where the mass loss rate 
is fixed to the SFR. 

6. Finally, the quantity az is the ratio of the metallic- 
ities of infalling gas and the ISM: Zi^jZisu- The 
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inclusion of az allows for enriched inflows (pos- 
sibly from previously ejected material now being 
re-accreted from the galaxy halo). Additionally, it 
appears in the model for SFR as a simple way to 
quantify the recycling of enriched halo gas into the 
ISM (i.e. an additional source of fuel for star for- 
mation). The parameterization of a^ is taken from 
momentum-d riven wind simulations in iDave et al.l 
(l201lE[20T2h : 



az = (0.5 - O.lz)(M,/lO^"M0) 



\0.25 



(9) 



(although they note that it is a "crude" parameter- 
ization) . While we aim to test models beyond those 
that are momentum-driven, a^ is not reported for 
these simulations. Therefore we use this parame- 
terization for all of the models with winds. For the 
no wind case we consider a^ = 0, since there can be 
no recycling of previously ejected material in this 
scenario. We consider the effects of modifying a^ 
in the discussion below. 

7.2. Comparison of equilibrium models to 
intermediate-redshift data 

Figure |8] shows how the equilibrium models com- 
pare to the intermediate-redshift MZ and SFR-M* re- 
lations. In order to get the best constraints on the 
models, we combine our sample with binned DEEP2 
data representing the MZ relation for z ^ 0.8 galax- 
ies (|Zahid et al l I2OII D and the SFR-M , relation for 
2; ~ 0.5 - 0.7 galaxies (jNoeske et al.l 120071 ). As we men- 
tioned above, since metallicity calibrations and nucle- 
osynthetic yields are uncertain, wc fit for the yield that 
reproduces the MZ data. The yields that minimize y^ 
(including uncertainties in both metallicity and masso) 
are y = 0.012,0.013,0.010 and 0.008 for the momentum 
driven, energy driven, constant wind (77 = 1) and no wind 
(77 = 0, ttz = 0) models respectively. We next describe 
how each model compares to the data: 

No wind — Not surprisingly, equilibrium models that ex- 
clude outflows are unable to reproduce the data. Most 
drastically, when 77 = and az = 0, Equation |4] has 
no mass dependence. Hence, the predicted MZ relation 
is implausibly flat, and the slope of SFR-M, relation is 
also too flat. Figure [5] shows that, without winds low 
mass galaxies make too many stars (and metals), and 
are unable to get rid of the metals that they make. It is 
worth noti ng that a^, is not s trictly zero in the no-wind 
simulation (jDave et al.ll2011a|) . since infalling gas can be 
enriched by other galaxies or could represent the accre- 
tion of a satellite galaxy. A non-zero Uz would further 
increase the SFR and amplify the disagreement between 
the SFR-M* data and the no-wind model. 

Constant Wind — A model where the mass-loss rate is 
equivalent to the SFR shows good agreement with the 
SFR-M, relation at all stellar masses. However, this 
model predicts an MZ relation which is too flat at low 
stellar masses. 

9 We take x? = {OH, - /{M,))V('^Sh,» + /'(A^O^'^i,,), where 
f{Mi) and f (Mi) are the model and its derivative evaluated at 
Mi . x^ is the sum over the i observations of x? • 



Momentum Driven Wind — Of the models shown in Fig- 
ure [51 the momentum driven wind shows the best match 
to the MZ data. However, the metallicities predicted by 
this model turn up at higher masses, owing to the de- 
pendence on Uz ■ Equation 2] shows that the predicted 
metallicities diverge as a^ approaches unity at higher 
masses. While the disagreement between the model and 
data is only apparent in the highest mass bin, the upturn 
remains impla usible. Locally, the M Z relation flattens to 
higher masses (jTremonti et al.|[2004[ ). so the model shown 
in Figure H] crosses the z ~ 0.1 relation at M ~ 10^°'^ 
Mq. This o bservation suggest s that the mass scaling for 
az given bv IDave et al.l ( 20121 ) may be too steep at high 
masses. 

In addition to the MZ relation, we compare the 
momentum-driven wind model to the SFR-M, relation 
in the top panel of Figure [51 In this case, the model 
slightly under-predicts the SFRs for the galaxies in the 
present study as well as for the lower mass DEEP2 mea- 
surements. This discrepancy between the data and mod- 
els at low masses is qualitatively simi lar to what i s found 
at higher mass es and l ower rcdshifts. IDave et"al] ()2011aD 
and Weinmann et al.l ([2012) find that both simulations 
and semi-analytic models are unable to reproduce the 
population of low mass galaxies with high specific SFRs. 

Given the uncertainty in Uz , we next consider whether 
modifying this variable can produce a momentum driven 
wind model that is a good fit to both the MZ and SFR- 
M» data. If we remove the mass-dependence of az in 
Equation [9l(az oc A4^°), the predicted SFR-M* relation 
appears very similar to the constant wind model (a good 
match to the data), but the predicted MZ relation (not 
shown) is far too fiat. In summary, modifying Uz does 
not improve the agreement between the data and the 
momentum driven wind model. 

Energy Driven Wind — The energy driven wind models 
shown in Figure [51 are too steep compared to both the 
MZ data and the SFR-M* data. However, as we al- 
ready noted, some of this steepening of the MZ rela- 
tion occurs because az may be too large at high masses. 
In fact, an energy driven wind model with az oc M^'^ 
and y — 0.022 appears very similar to the momentum 
driven wind model shown in both the top and bottom 
panels. For this model, the MZ relation model reduces 
to a straight line that runs through all of the data, and 
SFRs are slightly under-predicted at low masses. Ul- 
timately, we cannot distinguish the momentum driven 
model from an energy driven model with a^ oc M^'^. 
However, we do note that this variable shows some mass 
dependence in the momentum-driven wind simulations 
([Dave et al. 2011a,b), so it unlikely that this dependence 
is non-existent for energy-driven winds. We therefore 
conclude that energy-driven winds are unlikely to explain 
the data. 

Recent simulations from Hopkins et al. (2012) have 
suggested that outflows depend strongly on gas surface 
density, so that winds from low density dwarf galaxies 
are dominated by energy-driving from supernovae and 
stellar winds. More massive star-forming galaxies with 
higher gas densities, on the other hand, are predicted to 
prefer momentum driven winds. The present data do not 
indicate a mass-dependent 77; Figure [H shows no evidence 
for steeper MZ and SFR-M* slopes towards lower masses. 
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This analysis (and Figure |5]) shows that there may be 
a tension between the data and the equihbrium model. 
On one hand, models which best match the MZ re- 
lation under-predict the amount ol star formation in 
low mass galaxies (and even intermediate masses probed 
by DEEP2). Likewise, models that best reproduce the 
SFR-M* relation over-predict the metallicities of galax- 
ies. This trend can be understood as a generic prop- 
erty of the models, because both Equations [3] and |4] are 
proportional to (1 + ri)~^{l — az)~^; changes in these 
model components alter the MZ and SFR-M* relation in 
the same manner. However, if preventive feedback from 
gas heating has been overestimated (C underestimated in 
Equations [3] and [5]), then SFRs can b e increased without 
modifying the MZ relation. In fact, iDave et"an (|2012[ ) 
note that Cgrav is taken f rom simulations that do no t 
include metal-line cooling (|Faucher-Giguere et al.|[20TTI ). 
and raise the question of whether this effect is important. 
Alternatively, if ( is correct, then the data may suggest 
a deviation from the equilibrium model. 

Ultimately, we conclude that this disagreement re- 
mains subtle. While we have argued that emission line 
selection does not largely bias our results, we cannot 
rule out a small effect. As an example, if the SFRs pre- 
sented in Figure [5] are high by 0.1 dex, we would expect 
an unbiased sample to have an SFR-M* relation that 
agrees better with the momentum driven wind model at 
low masses. Additionally, according to the FMR, SFRs 
which are 0.1 dex lower should be accompanied by metal- 
licities which are around 0.03 dex higher, in even better 
agreement with the momentum driven wind MZ relation. 
Even without emission-line selection bias, it is important 
to recah that the [O II], H/3, and SED-based SFR indi- 
cators discussed in the Appendix differ systematically by 
up to 0.2 dex. However, modifying the SFR under these 
scenarios can not satisfactorily bring the data in line with 
the momentum-driven model. Even in the higher mass 
data alone, this model predi cts a steeper SFR-M ^, rela- 
tion than the DEEP2 result (jNoeske et al.ll2007[ ): there- 
fore, the discrepancy which we have identified is apparent 
in a magnitude-limited sample. We conclude that Ha fol- 
lowup spectroscopy of our sample (combined with a mass 
selected control sample) can definitively assess the effects 
of emission line selection, while simultaneously reducing 
errors on dust, SFR, and metallicity. We leave this work 
as the subject of a future paper. 

Finally, it is worth pointing out that systematic 
uncertainties in our metallicity calibration cannot ex- 
plain the discrep a ncy b etween the data and models. 
iKewlev fc Ehisonl pOOl) show that local (SDSS) MZ 
relations have similar lo w mass sl opes under most 
calibrations. The .Trem onti et al.l (|2004( ) calibration 
and the direct-method e lectron temperature calibration 
([Andrews fc Martinill2013i) are two exceptions- both pro- 
duce MZ relations that are steeper than what we have 
used here. Hence, if we calculated metallicities in the 
same way as Tremonti et al. or by using electron temper- 
atures, the tension between the data and models would 
be ampli fied. On the ot her hand, the calibrations re - 
ported bv lPilvuginl (pOOll ) and lPilvugin fc ThuanI (pOOl 
give MZ relations that are flatter than most. However, 
we note that they also have metallicities around 0.5-0.7 
dex lower than those from KK04. These lower metal- 
licities would introduce other difficulties, as they would 



require normalizin g yields that fall below the plausible 
range indicated in iFinlator fc Pavel (|2008D . 



8. CONCLUSIONS 

In this work we have placed the first constraints on the 
low mass, intermediate redshift MZ relation. By using 
emission line selection in the COSMOS field, we are able 
to efficiently identify and measure the metallicities from 
26 galaxies reaching masses of 10® M© at 2; '--^ 0.65. Com- 
bined with previous measurements from magnitude lim- 
ited samples, these data extend our knowledge of the MZ 
relation to masses an order of magnitude smaller. This 
limit is comparable to the low mass limit of the local MZ 
relation determined from SDSS data. Therefore, for the 
first time we are able to measure the metallicity evolution 
of intermediate redshift galaxies at M < 10^ M©. Com- 
pare d to the z ~ 0.1 relation reported bv lTremonti et aLl 
()2004i) . we find an average metallicity that is 0.12 dex 
lower for galaxies with 10®-^ < M/Mq < lO''". We show 
that this measurement is consistent with an MZ relation 
that evolves more strongly at low stellar masses. We in- 
terpret this mass-dependent evolution as consistent with 
downsizing trends, where higher mass objects have less 
leverage to alter their gas-phase metallicities after most 
of their stars have been assembled. 

An important development in our understanding of 
metallicity evolution is the discovery that the scatter in 
the relation can be reduced by accounting for star forma- 
tion. This measurement was quantified as a plane, and 
dubbed the Fuii d ameii tal Metallicity Relation (FMR) by 
iMannucci et al.l ()2010f ). Using the present sample, we 
find that a planar relation does indeed exist among low 
mass, intermediate redshift galaxies. Consistent with the 
findings of Mannucci et al., we see no evidence for evo- 
lution of the FMR, as our emission-line selected sam- 
ple falls in good agre ement with their z ~ 0.1 FMR 
(jMannucci et al.l 120111) . Furthermore, while the lowest 
mass galaxies in our sample show significant scatter with 
respect to the FMR, we attribute this scatter to mea- 
surement uncertainties for faint galaxies. Using a Monte 
Carlo simulation, we determine that the intrinsic scatter 
is consistent with zero, and should be less than 0.16 dex 
(95% confidence) at M < lO^-^ Mq. 

The MZ relation is an important probe of galaxy evolu- 
tion models, as accretion and galactic outfiows modulate 
the gas-phase metallicities of galaxies. Hence, compari- 
son to model predictions can help us to gain insight into 
galaxy formation. We have combined our MZ relation 
with the SFR-M* correlation that is also measured from 
our data, taking higher mass d ata from the literature 
("Noes ke et alll2007l: IZahid et al.ll201lD. We com pare to 
the family of models outlined in lDave et al.l ()2012l ) , where 
it is assumed that galaxies prefer to maintain an cqui- 
librium between inflow s, outfiows, and star formation 
(IFinlator fc Dav e 2008; Bouche et al.ll2010l:lDutton et al.l 
120101: IDave et al.ii2012i) . We find that models which pre- 
dict the MZ relation may under-predict the SFRs of low 
mass galaxies, and at the same time, models that pre- 
dict the SFRs of low mass galaxies tend to over-predict 
their metallicities. While this finding could represent a 
breakdown of the equilibrium model in low mass galax- 
ies, it could alternatively be an indication that feedback 
from gas-heating has been overestimated in simulations 
(i.e. .Fauchcr-Gigucrc et al., ,2011.) . To solidify this re- 
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suit, we have begun Ha and [N II] spectroscopy to fol- 
lowup the present sample and simultaneously measure a 
mass-selected control sample. These data will clarify the 
effects of emission line selection and greatly reduce the 
statistical uncertainties in our metallicity, dust, and SFR 
measurements . 

To conclude, these observations have provided a valu- 
able look at the metallicity evolution of low mass galaxies 
outside the local universe. We look forward to the physi- 
cal insights and new constraints that can be gained from 
larger samples and a better characterization of system- 
atic uncertainties. 
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APPENDIX 
ASSESSING SFR MEASUREMENTS 

In [jni and iJT] we draw conclusions from SFRs that we derived from SED fits; therefore, it is important that we 
use rehable SFRs. In order to assess systematic uncertainties, we compare this measurement to two additional SFR 
indicators: II/3 and [O II] . For the H/3 method, we use the dust exti nction from the SEP fits, scaled up by a factor of 2.3 
to account for the difference between stellar and nebular extinction (jCalzetti et al.l[2000l ) . We also correct the H/3 fluxes 
for the small amount of stellar absorption, as discussed above. Next, the H/3 luminosities are converted to Ha, assuming 
a B aimer decrement of 2.8. We then use the Ha-SFR calibration given by Kennicutt (1998), scaled appropriately to a 
iChabrieii p003) IMF. To derive SFRs from [O II] luminosity, we use the calibration given bv lMoustakas et al.l ()2006f) , 
which includes an Mb dependence, but does not require that the line measurement be corrected for dust. 

In order to make a direct comparison to the SED-derived SFRs, we must correct the emission line fluxes for both slit 
losses and and extraction aperture losses. Rather than model these losses (which would not account for uncertainties in 
the absolute flux calibration), we compare our data to narrowband imaging in C OSMOS. First, we i nfer the emission 
line flux from the narrowband (NB816) photometry, using the methods outlined in lLy et aD (|2007t ) and'Taka hashi et al.l 
([2007). Since the NB816 bandpass has a non-uniform throughput, we calculate a correction to this flux based on the 
relative throughput at the observed wavelengths of the emission lines. (This correction can be important, since it is not 
unusual that both of the [O III] doublet lines are found away from the central wavelength of the NB816 filter.) Because 
the COSMOS catalog contains total magnitudes (under the reasonable assumption that the galaxies are unresolved 
with 1.5" FWHM spatial resolution), the line fluxes inferred from narrowband imaging are not subject to aperture 
losses. For 17 galaxies in our sample with emission lines that are bright enough to be detected in the narrowband 
imaging, we can infer the aperture losses by comparing these narrowband imaging fluxes to our spectroscopic fluxes. We 
find typical correction factors of 1.1-2.0, with a median of 1.4. We adopt this average correction to explore systematic 
differences between star formation indicators. 

Figure [9] compares the SFRs derived from the SEDs, H/3, and [O II]. In addition to the SED fits described in §3.2 
(which use exponentially declining star formation histories with r as a free parameter; red points), we also make the 
comparison for constant star- forming models (blue points) . This exercise shows that the SFRs derived using different 
methods are correlated, so that we can easily distinguish the galaxies with high and low SFRs. However, Figure [S] 
also shows that there are systematic shifts as large as 0.2 dex between the different estimators. Determining the most 
accurate star formation indicator is beyond the scope of this paper. Ultimately, a systematic uncertainty of 0.2 dex 
does not affect our conclusions. In [J6l we show that the scatter in the MZ relation is correlated with SFRs. As with 
the local FMR, metallicity depends only weakly on the SFR, so that a 0.2 dex offset in SFR translates to a 0.03 dex 
shift in the metallicity predicted by the local FMR ([Mannucci et al. 2010, 2011). This difference in metallicity is much 
smaller than the uncertainties on the FMR residuals shown in Figure [T] Likewise, in S}7]we use our SFR-M, relation to 
constrain models. However, 0.2 dex of systematic uncertainty is small compared to the dynamic range of SFRs shown 
in Figure [8] An systematic offset is unlikely to alter our conclusions, as the SFR-M* trends discussed in ^can be 
seen (albeit less obviously) from the higher mass data alone. 
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Figure 9. Star formation rates derived tliree different metfiods are compared. Blue points assume a constant star formation history, and 
the red points assume exponentially declining SFRs (t models) discussed in §3.2 The H/3-derived SFRs are corrected for dust extinction 
using these SED fits. (The center and right panels— where the SED-independent [O II] SFRs are compared- show that the r models tend 
towards slightly less dust.) The [O II] and H/3 SFRs are corrected by an average factor of 1.4 to account for slit-losses. Systematic shifts 
between the different indicators are at most 0.2 dex. 



